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Abstract 



We study the flavorless gluon propagator in the Landau gauge from high 
statistics lattice calculations. Hypercubic artifacts are efficiently eliminated 
by taking the Y^P^i — ► limit. The propagator is fitted to the three-loops 
perturbative formula in an energy window ranging form ~ 2.5 GeV up to ~ 
5.5 GeV. a s is extracted from the best fit in a continuous set of renormali- 
sation schemes. The fits are very good, with a x 2 P er d.o.f smaller than 1. 
We propose a more stringent test of asymptotic scaling based on scheme in- 
dependence of the resulting A^. This method shows that asymptotic scaling 
at three loops is not reached by the gluon propagator although we use rather 
large energies. We are only able to obtain an effective flavorless three-loops 
estimate A^L = 353 ± 2+iq MeV. We argue that the real asymptotic value 
for A^jg should plausibly be smaller. 
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The non-perturbative calculation of the running coupling constant of QCD is certainly 
one very important problem. This program has been performed using the Schrodinger 
functional the heavy quark potential 0|§, the Wilson loop Q, the Polyakov loop || 
and the three gluon coupling An usual approach to determine the strong coupling 

constant (for instance in ref. 0) consists on the evaluation of an ultraviolet quantity and 
on the further comparison of the numerical data with perturbative predictions, the gluon 
propagator at large momenta is a good candidate to be used. 

Much work has been recently devoted to the study on the lattice of the gluon propagator 
the effort being mainly concentrated on infrared behavior. The authors of |TD[ have 



started an ultraviolet study by comparing the asymptotic behavior to the one loop QCD 
prediction. In this paper we would like to concentrate on the asymptotic behavior of the lat- 
tice propagator in the Landau gauge by a systematic use of the three-loops QCD prediction, 
leaving a study of the infrared to a later publication. 

The propagator is evaluated to a much better statistical accuracy than the coupling 
constant computed in and 0. The three-loops fit should allow a determination of the 
strong coupling constant and hence of Aqcd to a good accuracy in view of the very large 
statistics we have accumulated (1000 configurations at /3 = 6.2 on a 24 4 lattice). This can 
be done in any renormalisation scheme as long as the anomalous dimension of the gluon 
field and the beta function are known to three loop. A test of scheme dependence can thus 
be performed rather extensively. It will turn out that although the statistical accuracy is 
very good as expected, there remains a large systematic error due to the non asymptoticity 
of the gluon propagator. 

In section 1, the general principle of the method is explained. In section 2 we perform 
a discussion of the lattice artifacts related to the hypercubic geometry, and we propose a 
solution to eliminate them. In section 3 the fit is performed in the MOM scheme. In section 
4 we discuss in general the scheme dependence. We conclude in section 5. 

I. GENERAL DESCRIPTION OF THE METHOD 

The Euclidean two point Green function in momentum space writes in the Landau gauge: 

G?ir 2 fe -p) = G^\v 2 )5 aia2 (<w 2 - (i) 

where a±, ci2 are the color indices ranging from 1 to 8. 

In any regularization scheme (lattice, dimensional regularization, etc.) with a cut-off A 
(a" 1 , (d — 4) _1 ) the bare gluon propagator in the Landau gauge is such that 

lim dln^GgL^A)] 
A^oo dlnp 2 

is independent of the regularization scheme^} Lattice calculations provide us with an evalu- 
ation of the bare propagator, and hence of the logarithmic derivative in eq. (|2|). The latter 



3 We will return to this statement in the next subsection. 
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is an observable on which we will concentrate in this paper and from which we will compute 
the strong coupling constant. 

In the MOM (or MOM) scheme, calling Z 3 (q) the gluon renormalisation constant: 

Z 3 (q,A) = q 2 G?l e (q 2 ,A) (3) 

where q is a positive energy scale defined by q 2 = p 2 , it happens that the expression in eq. 
(|) is simply equal to liniA-^oo din Z 3 (q, A)/d\nq 2 . 

An important point has to be stated here. The fact that Z 3 is the renormalisation 
constant in the MOM scheme does not constrain us to stick to the MOM scheme all through. 
Equation (0) defines unambiguously a quantity, the evolution of which we can study in any 
scheme. This evolution is given by the anomalous dimension of Z 3 . 

dlnZ 3 (q,A) f 7o 71 2 72 3 . 4 A . . 

_ — . — = r(ct) = - —a + — —a + 7— -r^« + 0[a ) (4) 

ding 2 \4:7i (47r) 2 (47r) 3 J 

where it is understood that the coupling constant in a given scheme is a function of q such 
that 

& a of \ A) 2 Al 3 As 4 , 5n /r\ 

— — = plot) = or — - — —a — - — —a + Ola ) 5 

<91ng Hy ! 2tt (2tt) 2 (4vr) 3 V ; V ; 

with 

13 

A, = 11, /3 a = 51, 70 = — (6) 

in the flavorless case, 71,72 and (3 2 being scheme dependent. As we shall see later, there 
is one scheme-independent relation between 71,72 and fa- To be specific, in the flavorless 
MOM scheme, for which 02 is known [I7J : 

29 

/3 2 ^ 4824., 71 = —, 72^ 960. (7) 

o 

The calculation of 71 and 72 will be explained later. 

From eqs. and (|5]) it is easy to integrate simultaneously, up to three loops, In Z 3 (q) 
and a(q) provided one is given the values Z 3 (fj,) and at some initial value q — 



q(a) = [i 



a 32ty 2 (3 + 16iT(3ia(jj) + (32a 2 (/j,) 



a{jj) v 327T 2 /? + 167r/5ia + /3 2 a 2 



e W {^(--^) +7 ^((3 (3 2 -4p 2 )[H(a)-H(a( f ,))}\ , (S) 



and 



Z 3 (a) = Z 3 (fi) 



a 



2H , 72 70 

"0 / 32tt 2 /? + 167r/3ia + (3 2 a 2 \ * 2/3 ° 



^a(fi) J \327r 2 (3 + 16irf3ia(n) + f3 2 a 2 (jj,) 

A) A 



exp { ( 271 - ^ - 4^ [tf(a) - if («(/,))]} ; (9) 



where, for 2/3 /3 2 - 4/3? > 0, 



H(a) = —j^^=^^= arctan ( = J^ L- I . (io) 

/2/3 /3 2 - 4^? V 47r V / 2AA I 4^ 



Eqs. (HH) give a parametric representation of the exact solution of the coupled differential 
equations ([|) and (0), where a can be considered just as a parameter connecting Z 3 and g. 
The use of this parametric representation allows to work with exact solutions of Eqs. @) and 
(H), the computation of a(q) from Eq. ([8]) being only approximatively possible. However, 
as explained in section 4, no important discrepancy comes from the consistent perturbative 
expansion on a of the exact solutions (|8],|9[). 

Our general "strategy" will be to fit the lattice results with one solution ([|||). If an 
acceptable best fit exists, it will prove that the lattice results do agree in the considered 
energy range with three-loops perturbative QCD. Moreover, it will provide us with the 
best initial values Z{p) and a(//). Z(p) is just an overall multiplicative constant, but the 
knowledge of a(/x) in a given scheme allows to compute an estimate of Aqcd in this scheme, 
and hence Aj^g. 

To our knowledge, this method of determining the strong coupling constant on the lattice 
is new. In |10| , the asymptotic behaviour of the gluon propagator has been compared to the 
one loop perturbative QCD prediction. The authors indeed find rough agreement. However, 
the value of A that one can extract from the one loop fit is not meaningful since, 70 being 
scheme independent, one does not know in which scheme it is computed. In the following we 
will systematically compare the two loops to the three loops result in a large set of schemes, 
and we shall conclude that at the accessible scales of 2.5 - 5. GeV, the gluon propagator 
does not scale to two loops and hardly does to three loops in the most favorable schemes. 
In the MS scheme, no sign of scaling is found. 

A. Three loops expansion 

In a general renormalisation scheme, which we call SC, and in a fixed gauge 

41(P 2 ,A) = Z 3SC (q,A)[G^(p 2 , q ) + 0(1/ A)] (11) 

(2) 

where Z 3 gc is the renormalisation constant in the scheme and G^q the renormalized gluon 



propagator. Notice by the way, that from eq. (|TT|) one immediately obtains the well known 
above-mentioned result that 



Um d\nZ 3SC (q,A) = d\nG^(p 2 ,q) 
A^oo dlnq 2 ding 2 



4 We write the formal solution of Eqs. (§J5|) for the case of positive discriminant because this is 
the case, for instance, of MOM and MS schemes. 
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is finite and independent of the regularization scheme, while the r.h.s is independent of p 2 . 
From (|) and (|rfj) 

Z 8 (g, A) = ? 2 41(<? 2 ; A) = Z 3SC (<7, A)[g 2 4?(? 2 ,?) + 0(1/A)] (13) 

whence 



Hm d\nZ 3 (q,A) = Um d\nZ 3SC (q,A) + d\n[q 2 G^(q 2 , q)] , { { ) 

A^oo cilng 2 A^oo dlnq 2 dlnq 2 



Expanding eq. (|14D in a introduces a dependence on the scheme in which a is expressed. 
We will now specify the scheme SC to be the MS scheme since the gluon propagator anoma- 
lous dimension has been computed to three loops (see eq. (8) in [|T5|j): 

.. d\nZ 3m (q,A) 13 531 2 29311. 3 4 , 

hm . ' ~ ; — -a ; — rrra 2 ; — --a 3 + 0(a 4 ) (15) 

A^oo ding 2 2(4vr) 8(4vr) 2 32(4tt) 3 v ; v ; 

where a is the coupling constant in the MS scheme and where the e = d — 4 — > limit has 
been taken. Using now0 the expression for q 2 G^(q 2 , q) in |16| we can rewrite eq. ( |14|) as 

d\nZ 3 (q,A) _ 70 _ 7 X _ 2 % _ 3 4 

hm — - — = T(a) = -j—^-a - - , cT - - ,, cT + U{or) 

A^oo ding 2 v y (4tt) (47r) 2 (4tt) 3 v 7 

13 155.3 9 6656. n Ax ,„„,. 

-a - - — -a 2 - - — -a 3 + 0(a 4 ) (16) 



2(4vr) (4vr) 2 (4vr) 



B. General three-loops schemes 

If the strong coupling constant is expressed in some other scheme, which we shall call 
generically the "tilde" scheme, once known a as a function of a, changing scheme simply 
amounts to a change of variables. If 

a = a + —a + — — --a + ... (17) 

47T (47r) 2 



then as shown for example in (17 



A m = Ae p 2 = p 2m + 2[3 (b - a 2 ) - Aafc . (18) 
The 7i's being defined from 

dlnZ 3 (g,A) _^p_ _ _jYi_ 2 _ 3 Q{ 4) ( ) 

A™ ding 2 W ~ (4vr) (4vr) 2 (4tt) 3 + K h { ' 



5 [q 2 G^(q 2 ,q)] 1 is equal to J ren as given in eq. (8.13) in [16 |. 
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are then computed from the change of variables (|T7|) inserted into (|16|), finally leading to 
the relations 



A exp 



7i ~ 7i 
27o/5o 



(20) 



02 — Ams + 2/^o 




-4A 



7i ~ 7i 
7o 



(21) 



It results that fc, 7i, and 72 are not independent parameters. Knowing two among these 
three parameters allows the computation of the third one. We stress that in any "tilde" 
scheme, only two parameters are required for any purpose we need in the study that we 
present in this paper. In other words, the set of renormalisation schemes is a two parameter 
space in which the MS is one point, and MOM another. On the following, we will use 
the coordinates (71, 72) to characterize any renormalisation scheme of the two parameter 
space. The possible use of this large parameter space of schemes will provide us with a 
tool which we will exploit extensively in the following. Note that when 71 = 72 = 0, the 
coupling constant a is nothing but the "effective charge", flSfl, associated to the observable 



defined in eq. (Q). In reference [18] this charge was argued to be the proper object on which 
perturbation theory applies. We will generalise the latter philosophy by considering also 
schemes with a non- vanishing 71,72- 

In any "tilde" scheme, once evaluated 5(/x), we can compute A up to three loops. For 
a discussion of the different formulae, see HI7| . At two loops we will use the conventional 
formula 



A(c) 



/iexp 




Pi 

"0 1 

-'0 



(22) 



At three loops, we will use, when A = 2p p 2 — 4/3 2 > 0, either the unexpanded formula 



§1 



or the expanded one: 



arctan 



2(3 1 + (3 2 a/4ir / 



— arctan 



2/3i 



(23) 



16tt 2 /^ 



-a 



(24) 



In eqs (p3| - p4|) we have omitted for simplicity to write the /1 2 dependence of a. 



6 



II. HYPERCUBIC ARTIFACTS AND OTHER 0{a 2 p 2 ) EFFECTS 



We refer to @ for a description of the lattice simulations which have been performed, 
the calculation of the Green functions, their Fourier transform, the checks of the S aim color 
dependence of the propagators, and the set of momenta considered for the different lattices 
studied. We will use in this study 1000 configurations at (3 = 6.2 on a 24 4 lattice and 100 
configurations at (3 = 6.4 on a 32 4 lattice. The large statistics involved will reduce the 
statistical error to a negligible value. 




FIG. 1. Plot (a) shows (open circles) Z^{q), eq. ^), as a function of the scale q with a distinct 
point for each orbit. It also shows (black circles) Z^{q) with the factor q 2 replaced by q 2 in eq. ^). 
Even with the latter which reduces somehow the dispersion, the difference between the individual 
orbits at a same q exceeds by far the statistical errors. Plot (b) compares the "democratic" selection 
among orbits (lower curve), as a function ofq and the curve extrapolated to = as in eq. $2{\). 
The latter procedure gives obviously a much smoother result. 

In a finite hypercubic volume the momenta are the discrete sets of vectors 

= ( 25 ) 

where are integer and L is the lattice size. As explained in 0, we have averaged the 
propagators on the hypercubic isometry group H^. The momenta corresponding to e.g 
n M = (2,0,0,0) and = (1,1,1,1), belong to different orbits although they both have 
the same p 2 = 4(2n) 2 /L 2 , i.e. they belong to the same orbit of the continuum isometry 
group 5*0(4). Our statistical errors are so tiny that the difference between the evaluated 
propagators for two such orbits of same p 2 are quite visible, as shown in fig. |l](a). Such 
differences are understood as an 0(a 2 p 2 ) artifact of the hypercubic geometry of the lattice. 
For example, if one uses, rather than (|23|), the momentum 
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^ = -sin(^). (26) 



a V 2 

the resulting momentum squared differs by a relative 0(a 2 p 2 ): 

f = p 2 - -UV 41 + . . . , where p [A] = £>J (27) 

To reduce the hypercubic artifact, the authors of [10] advocate both the use of (^Bj) rather 
than (|25|), and the restriction to a "democratic" set of momenta, i.e. those for which p 2 is 
rather equally distributed on the different directions (pV). 




FIG. 2. Plot (a) shows (full line) the best fit of the three loops formula to G2(q 2 )q 2 with the 
"democratic" selection of orbits for 2.74 GeV < q < 4- 21 GeV. The fit is prolongated outside the 
energy window in dashed line. The x 2 P er d.o.f is 15.5. Plot (b) shows the best fit of G2(q 2 )q 2 
extrapolated to pM = according to eq. §2^) for 2.78 GeV < q < 4.32 GeV. The x 2 per d.o.f 
is 0.52. The improvement of the x 2 P er d.o.f by a factor ~ 30 is dramatic. Notice furthermore 
that the number of points is larger in plot (b) because the "democratic" selection eliminates all the 
orbits for some values of q. 

We will now propose another approach which we will eventually compare to the "demo- 
cratic" one. In general, different orbits with the same p 2 have different pM. For example 
(2,0,0,0) has nM = J>J = 16 while (1,1,1,1) has n 4 = 4. On the lattice, the function G^ 2 \p 2 ) 
defined in ([!]) is indeed a scalar form factor invariant under H4 that we shall assume to be 
a smooth function. The general structure of polynomials invariant under a finite group is 
known from invariant theory |19 |. For our purpose it is sufficient to know that any smooth 



if 4 -invariant function is indeed a function of the 4 invariants p™ = X^P^ ; n — 2,4,6,8. 
We will neglect the invariants with degree higher than 4 since they vanish at least as a 4 
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and parametrize the lattice two-point scalar form factor as a function (p 2 ,pM). When 
several orbits exist for one p 2 , it is possible to extrapolate to pM = o and we define 

G<l(p 2 ) ee hm G£V,P W ) (28) 

It is easy to see that, neglecting 0(a 4 ), if one uses ( |26"D instead of ([25"D, defining accordingly 
pM anc i (5|^(p^,p[ 4 l), one has 

hm G£ t V,pW) = hm Gf2(p 2 ,pt 4 l) (29) 

pW-f-0 pI 4 l-*0 

where the limit in the r.h.s. is taken at constant p 2 . We have indeed numerically checked 
form eq. fl2"9"|) the absence of sizable 0(a 4 ) effects. 

In fig. |l](b), the "democratic" Z 3 (q) and the one computed from eq. (|28|) are compared. 
The latter provides a much smoother Z%(q). This is confirmed by the x 2 of the fit in the 
next section. The best fit of the "democratic" Z3 using the eqs. @) and (||) gives a \ 2 P er 
d.o.f of 15.5, see fig 0(a), while the best fit to Z 3 computed from eq. fl2"gp gives a x 2 per 
d.o.f of ~ 0.52, fig. 0(b). For sure, the pM — > extrapolation increases the errors on Z 3 as 
compared to the errors in the individual orbits, but not too much. We have computed the 
errors on the extrapolated points using the jackknife method. Typically the extrapolation 
increases the errors by a factor ~ 2 which cannot account for an improvement by a factor 



~ 30 on the x P er d.o.f. While this paper was in writing appeared a study by J. P. Ma, flU" 
who suggests the use of the variable p 2 + a 2 p' 4 l /12 to cure hypercubic artifacts. He shows a 
significant smoothing (fig. 2 in ref 



We therefore conclude that eq. (28) allows to eliminate in a consistent way the hypercubic 
artifacts. This does not mean that all 0(a 2 ) artifacts are thus eliminated: for instance, the 
lattice artifacts oc a 2 (p 2 ) 2 , which do not break 5*0(4) invariance, are obviously not. Only a 
comparison of our results for different lattice spacings allows to estimate the latter effects. 



III. GENERAL FIT IN THE MOM SCHEME 

We first show the fit of our data at f3 = 6.2 (1000 configurations) with solutions of the 
coupled differential equations (^) and (|5p in the MOM scheme. It has been performed on 
the lattice propagator in the energy window^ 2.97 - 4.32 GeV. The result of the fit is 

Z 3 (4.32 GeV) = 1.625(5), 0^(4.32 GeV) = 0.3005(15), xVd.o.f. = 0.44 (30) 

The x 2 P er d.o.f. is significantly smaller than 1. This feature may be a sign of some 



correlation between the points at different values of the energy q. Using (p3j|2~2|) and ( |20| ) 
we obtain 

A§L ~ 0.346 A { ^ = 354.2 ± 2.5 MeV, A^l ~ 0.346 A { ^ ~ 454 MeV, (31) 

MS MOM ' MS MOM ' v ' 



-"This best fit has been performed on a slightly different energy window than in fig. 0(b). 
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where the error is only statistical. The difference between the two-loops and the three-loops 
result is large. We will discuss this feature in the next section. We have checked that the 
result is not significantly changed when other energy windows are used. 

The study of our data at (3 — 6.4 (100 configurations) will be useful to estimate lattice 
artifacts oc a 2 (p 2 ) 2 . The best fit is now obtained for the window 2.97 - 5.48 GeV . The 
results for such a fit are the following 

Z 3 (5.48 GeV) = 1.452(18), a^(5.48 GeV) = 0.255(4), xVd.o.f. = 1.1 (32) 

For comparison with (|30D, performing the fit in the energy window 2.97 - 4.32 GeV, we 
obtain a IV[ ^ [ (4.32 GeV) = 0.287(10). Notice that the _> o extrapolation meant to 
reduce hypercubic artifacts, eq. (|28|) , is slightly less efficient than at (3 = 6.2, maybe due 
to the smaller statistics. Still it allows a gain of a factor 15 for the x 2 with regard to the 
"democratic" analysis. It results 

AgL ~ 0.346 A (3 L = 346 ± 8 MeV, A^ ~ 0.346 A {c ^ ~ 441 MeV. (33) 

MS MOM ' MS MOM V ' 

The comparison between the results at (3 = 6.4 and at (3 = 6.2 seems to indicate that the 
lattice artifacts other than those oc a 2 p' 4 ' are not too important. We did not use our data 
at (3 = 6.0 since they do not reach a large enough energy scale. 



IV. SCHEME DEPENDENCE 

A. MS scheme 

It might look a little too involved to fit the data as a function of ol^x,(h) and then to 



(3) (3) 

convert A ~ into A^. Would it not be simpler to work all the way in MS, i.e. to perform 

MOM MS v 

the fit as a function of a(/i) and to convert directly the result into A^ ? We have done this 
exercise and found 

Z 3 (4.3GeV) = 1.623(5), a(4.3GeV) = 0.207(2), xV d -°- f - = °- 64 ( 34 ) 
leading via eq. ( p3|) to 

A® = 549 ± 9 MeV, A§1 ~ 765 MeV (35) 

IVI o IVI o 

The discrepancy between (^) and ([H]) comes as a big surprise, since we have only 
performed a change of variables. The only difference between the two methods comes from 
the truncation of the perturbative series. To be more precise, when we work with ol^^ to 
three loops, we truncate the series (|4]) and ([5]) beyond the third term. Changing variables 
to a means expanding a in terms of up to the third term, (|17| ), implementing this 

change into the functions r(a) and (3(a) and expanding the result up to the third term. 
This final expansion rejects in a consistent way some parts which are formally equivalent 
to four and higher loops. Could it be that four-loops effects are so large as to induce the 
discrepancy between (|35D and (31]), although we are working at an energy scale larger than 
4 GeV ? 
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FIG. 3. Plot (a) shows in full line 2F(a) in the MOM scheme as a function of Aj^, using 

eqs and (ffiJp. It shows in dashed line 2T(a) in the MS scheme as a function of A|^, using 

eqs and (p3[). The horizontal line corresponds to a typical value of the logarithmic derivative 

d In Zs(q)/d In q evaluated from our lattice data around 4 GeV. Plot (b) adds to the preceding plot 

the curves 2T(a) (resp. 2T(a)) truncated at first order i.e. 71 = 72 = (resp. 71 =72 = and at 

second order i.e. 72 = (resp. 72 = 0) in full lines (resp. dashed lines). In both cases the curves 

range from right to left with increasing number of terms. In all cases an energy scale of 4-1 GeV 

(3) 

has been used to compute A™.. 



To dig into this question, we have performed the following exercise. In fig. 0(a) we 
plot 2T(a) as a function of A^, using eqs (§]) and ([231). We do this in both schemes. The 
horizontal line, at -0.35 is about the slope din Z 3 (q)/d\nq evaluated from our lattice data 
around 4. GeV. What happens is now obvious, the curve corresponding to MS at three 
loops is far from the MOM one at such a large absolute value of the logarithmic slope as 

0. 35. It is also obvious from the fig. 0(a) that, had the slope been smaller in absolute value, 

1. e. had the horizontal line been, say at -0.20, the two schemes would have agreed much 
better. But the lattice results show that such a slope can only be reached at a much larger 
energy than 4. GeV even though we are working at three loops !! In other words, the lesson 
is that the Landau gauge gluon propagator reaches asymptotic scaling only at very large 
energies. We are indeed starting a study at larger energies. 

The next question to address is whether it is possible to choose the best scheme, in the 
sense of the scheme which is closest to asymptotics ? Of course we have a prejudice, from 
other theoretical and experimental sources about Aj^g, including our own analysis of a, |7|,|r7 



that the result in ( pTD is better. However, it would be more consistent to conclude about 
this question with the only evidences coming from the analysis of the gluon propagator. In 
fig. 0(b) we plot the same curve as in fig. 0(a), together with the equivalent plots, but 
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taking 72 = 0, and the ones with 71 = 72 = 0[|. It is clear that the three curves are much 
closer in the MOM scheme than in the MS one. This is an argument in favor of MOM while 
MS proves really bad for this problem: the two-loops is more than 200 MeV away from the 
three-loops, not to speak of the one-loop which does not even cross the line in our plot ! 



B. Search in the scheme space 

As long as we are engaged in comparing the schemes, why not exploit the richness of the 
large scheme space to which we can access up to three loops simply by defining arbitrary 
couple of parameters (71,72). From ([H]) one also knows 02 and from ( p0|) the ratio A^jg/A. 
Hoping to find schemes for which asymptotic scaling is reached, our strategy will be to scan 
the latter parameter space, and to try a fit of the data at — 6.2 in order to find the true 
asymptotic A^g. We select the domain of the "good schemes" by imposing some constraints 
which will be discussed right now, and the variation of A® in the latter domain will provide 
an estimate of one source of systematic uncertainty. 

In order to determine the good schemes we will impose three criteria. The first two are 
the following: 

l A H - A Sl < A — , SA < 5 max (36) 
with computed to three loops, using (£3), and A^ computed to two loops, i.e. taking 

,(3) A (3) I a (3) 



the same 71, but with 72 = 0, and (E2j); 5A is the difference |A^4 — A^l, with A^ 

' ~ 1 Mb exp Mb 1 7 Mb exp 



computed to three loops, using Q24|) . Typically A max will be taken around 100 MeV and <5 max 
around 10 to 25 MeV. The reason for these selection criteria should be rather clear: we wish 
the series of the functions /3(a) and T(5) to look like being in their convergence domain; the 
difference between the two-loops and the three-loops should thus not be exceedingly large 
as discussed in the preceding subsection. Finally we wish that expanding the formula for 
A^ 3 -* does not change drastically the result in order to keep some faith in the perturbative 
character of the latter formula. 

The criteria in ([3(]) are not restrictive enough to forbid exceedingly large values of 71 
and 72. To cure this we try some reasonable guesses about higher order terms. We extend 
to four loops the method which lead us to (|?p, 

3 „ , 2/32(71-71) 401(7! -7i) 2 (71 - 71 ) 3 

P3 — P3 H 1 9 3 

7o 7o 7o 

80o (71 ~ 7i)( 2 7i ~ 2 7i7i ~ 7o72 + 7072) 
7o 



7 We have kept for these curves the three-loops computation of A according to (|23| ) with non- 
vanishing 72: these are not genuine two or one- loop calculations. Later we will rather compare 
the three-loop result with the consistent two-loop calculation of the same scheme, using (|2^). Our 
qualitative conclusion about the "best" schemes will not change. 
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FIG. 4. Domain of the "good schemes" in the 71,72 plane. The level curves correspond to 
values 0/ Aj^g , ranging from ~ 380 MeV for the darkest down to ~ 345 MeV. The black circle 
indicates the position of the MOM scheme. The MS scheme is far outside. 



Wo 
'To 3 



(-57i + 67171 - 7171 + 5707172 - 3707172 - 2707172) 

_ 4 A) (73 -73) 
7o 



(37) 



where 7^ and /3j are in the MOM scheme, eq. (0), and where the beta function is defined up 
to four loops as 

of \ 00 2 01 3 02 4 03 5 /on\ 

Pick) = a — - — —a — -, — r^a — -, — 77 a . 38 

HK ' 2tt (2tt) 2 (4tt) 3 (4tt) 4 v ; 

Since 73 is unknown, we restrict the "tilde" scheme by imposing 73 = 73. We do not know 
/?3 in the MOM scheme, but we believe that it cannot be too large because the good fit 



of cv, extracted from the three gluon vertex in [|I7]], using the three-loops expression (i.e. 
with p*3 = 0), persisted down to rather low energies. Therefore, by bounding the distance 
between j3 3 and j3 3 we will keep /3 3 reasonably small. Now we are ready to add to P3) a 



third criterion to bound a reasonable domain of the "good schemes" : 

|&-/? 3 | <4tt/? 2 (39) 

It is worth noticing that this last criterion eliminates the larger values of accepted by 
the criterion (|36|). 

The three conditions ( j3~6| ) and ( |3~9"D define in the (71, ^-parameter space a domain of 
"good schemes" which is plotted in Fig. |j. The values of Aj^g vary between ~ 345 MeV and 
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~ 380 MeV. We will take this as an indication about our systematic uncertainty. However, 
this is not the whole story. All our analysis has shown that the gluon propagator is not 
asymptotic at three loops in the present energy range. The range ~ 345 MeV to ~ 380 MeV 
may be far from the real asymptotic value of but it provides an effective three-loop 

(3) 

A^J- . It is difficult to estimate the distance of the latter from the former, ignorant as we 

MS ' to 

are of the four- loops coefficients fa and 73. This distance can be large as we shall illustrate 
now. We plot in fig. |5](a) the behaviour of Aj^g, computed by obtaining the best fit over the 
strong coupling constant in the MOM scheme, as a function of the unknown fa assuming 
73 = 0. As can be seen, a positive^ fa would reduce the value of the four-loop effective 
Aj^g compared to the three-loop one. In fig. ||(b) we plot assuming fa = the behaviour 
of different schemes among our "good schemes" as a function of 73: 73 is computed from 
73 and fa = via eq. (|37|). We see that their predictions for Aj^g tend to converge for 
positive 73 where the value of A^g decreases. Of course convergence between the schemes 
is an indication that we are closer to the correct value of 73. 

Consequently, we expect that the asymptotic A 
which has been estimated in this work. 



(3) 

MS will be smaller than the effective A^ 



P=6.2 ; Volume=24 



.360 - 



.340 



.320 



.300 - 



.280 



' 1 I 1 1 1 I 1 1 
window=[2.97,4.32] (GeV) 
(a) 




0.0 



1.5 3.0 4.5 6.0 



A 
.420 

.370 

.320 

.270 

.220 

.170 



P=6.2 ; Volume=24 



window=[2.97,4.32] (GeV) 



0.0 



1.0 



2.0 



3.0 



FIG. 5. Plot (a) shows A^g as a function of the unknown fa in the MOM scheme assuming 
73 = 0. Plot (b) shows three different schemes as a function 0/73 assuming fa = 0. The solid line 
if for the MOM scheme. A tendency to converge is seen as 73 increases 



^The positive sign for fa is plausible as this is the sign of = 0..2 and as P 3 y^ > 0. 
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V. CONCLUSIONS 



Our initial aim has not been fulfilled. We cannot at the present energy scale give a reliable 
estimate of from the gluon propagator because it has not yet reached its asymptotic 
regime at three loops ! This is surprising since a scale as large as 5 GeV is often assumed, 
without further scrutiny, to be large enough to be asymptotic even at two loops ! The MOM 
scheme is among our "good schemes". It is amongst the schemes closest to asymptotic 
convergence. On the contrary MS is very far from convergence. 

We have learned an important lesson about the criteria of asymptoticity. It is not 
enough to have a good fit over a large energy range with an asymptotic formula to check 
asymptoticity. Our fits are good because the error we introduce by ignoring higher orders 
is only logarithmic. Our energy range, as large as it may look, corresponds to an increase of 
logg by only ~ 0.4. Therefore the higher loops effects can be mimicked by a simple rescaling 
of Aj^g. The difference in the functional behavior of higher orders would only be apparent 
on an energy range containing several "e-foldings" , i.e. changes of logg by several units. 
Such a study over several "e-foldings" up to a very high scale has been performed in 

Combining our results for /3 = 6.2 and (3 = 6.4, we end up with the following value for 
the effective three-loops estimate: 

A|| = 353 ± 2±?g MeV (40) 

where the first error is statistical and the second is the systematic uncertainty estimated 
from the scheme dependence. 

From a study of the plausible effect of the fourth loop, we have argued that the real 
asymptotic A^ should lay below the result in eq. fllPf) . 

We have discovered that the scheme independence of the result is a much more demanding 
criterion of asymptoticity than the quality of fit. We have used many "ad hoc" schemes 
simply defined by a couple of parameters (71,72)- This technique could be and should 
be extended to other methods to compute a s and more generally to any program which 
performs a matching to perturbative QCD, as long as several renormalisation schemes can 
be used. 

Finally we have proposed a new method to eliminate the hypercubic lattice artifacts, 
namely to take the limit pM — > 0, eq. (f^J). We demonstrated that this method is very 
efficient. 

From the result of the present analysis we have undertaken a lattice calculation at smaller 
lattice spacings, hoping to reach the asymptotic regime. 
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